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Abstract: Starting with two coupled Bethe-Salpeter equations for the quark-anti- 
quark, and for the quark-glue-antiquark component of the quarkonium, we solve the 
bound state equations perturbatively. The resulting admixture of glue can be partially 
understood in a semiclassical way, one has, however, to take care of the different use 
of time ordered versus retarded Green functions. Subtle questions concerning the pre- 
cise definition of the equal time wave function arise, because the wave function for the 
Coulomb gluon is discontinuous with respect to the relative time of the gluon. A striking 
feature is that a one loop non abelian graph contributes to the same order as tree graphs, 
because the couplings of transverse gluons in the tree graphs are suppressed in the non 
relativistic bound state, while the higher order loop graph can couple to quarks via non 
suppressed Coulomb gluons. We also calculate the amplitude for quark and antiquark 
at zero distance in the quark-glue-antiquark component of the P-state. This quantity is 
of importance for annihilation decays of P-states. It shows a remarkable compensation 
between the tree graph and the non abelian loop graph contribution. An extension of 
our results to include non perturbative effects is possible. 
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1 Introduction 



It is generally agreed that a quarkonium state is not only made up of a quark and an 
antiquark (QQ), but contains also a component with glue (QAQ). There are, in addition, 
other contributions like admixtures of several gluons or light quark- antiquark pairs which 
are not considered here. We emphasize that we are not talking about hybrids, where the 
glue component is an essential part of the wave function, but about the admixture of 
glue in ordinary quarkonia. We are mainly interested in equal time wave functions here, 
not, say, in light cone wave functions. 

The presence of a component with glue in the P-states has become of great impor- 
tance due to the work of Bodwin, Braaten, and Lepage pQ. Although this component is 
suppressed, it allows the QQ pair in the QAQ component to be in a color octet S'-state. 
The 5- state annihilation of this pair is of the same order as the QQ annihilation in 
the leading QQ component of the P-state, because the latter is also suppressed by the 
vanishing of the wave function at the origin. Logarithmic divergences cancel in both 
contributions. Similar considerations hold for P-wave production. 

In spite of the principal as well as phenomenological importance of the QAQ admix- 
ture it is hard to find quantitative statements about this part of the wave function in 
the literature. We will approach the problem in a systematic way here. Even if one is 
interested in the equal time wave functions only, it is necessary to start with the general 
wave functions involving different times. The reason is that quark and gluon propagators 
also depend on time differences, therefore a treatment involving only equal time wave 
functions appears not possible within our approach. We leave open the question whether 
our results could also be obtained within an equal time Hamilton approach. In any case 
some of the subtleties mentioned below can only be understood if one knows the general 
wave function with different time arguments. 

In sect. 2 we use as starting point two Green functions, one involving QQ and the 
other one QAQ. In the spirit of the usual derivation of the Bethe-Salpeter equation, but 
with more necessary technicalities, we define several irreducible kernels. These allow a 
decomposition of the Green functions and a derivation of two coupled Bethe-Salpeter 
equations for the QQ wave function and the QAQ wave function. 

In sect. 3 we proceed by using the framework of non relativistic QCD (NRQCD) 
[2] , which is appropriate for the case that the velocities of the quarks are small (order 
a s <C 1). We furthermore choose the radiation gauge VA = which is known to be 
the gauge best suited for treating bound state systems. In leading order of perturbation 
theory the QQ kernel becomes static and the QAQ wave function does not enter into 
the equation for the QQ wave function. The latter can thus be solved in the usual way. 
After this has been done one can insert the solution for the QQ wave function into the 
equation for the QAQ wave function and calculate the latter in lowest order. 

The result shows some delicate features. The wave function with the Coulomb gluon 
is discontinuous at t = 0, where t is the relative time of the gluon with respect to 
the time of the QQ system. It vanishes for positive t, while for negative t it can be 
simply understood as the Coulomb field moving along with the quarks. The tree graph 
contribution to the transverse gluon differs by a factor 1/2 from what one would expect 
from classical electrodynamics. The origin of this discrepancy can be understood, it is 
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due to the use of time ordered versus retarded propagators in the different approaches. 
For the transverse gluon there is a non abelian loop contribution, in which the transverse 
gluon couples to two Coulomb gluons while the latter couple to the quark and antiquark. 
This loop graph contributes to the same order as the direct coupling of the transverse 
gluon to quark or antiquark, because the direct coupling is suppressed by a factor of the 
order of |p|/m ~ a s . 

In sect. 4 we finally calculate the amplitude for QQ at zero distance in the QAQ 
component of the P-state. As stated above, this quantity is of importance for annihila- 
tion decays and production of P-states. The result cannot simply be written in terms 
of, say, the derivative of the QQ radial wave function at the origin, it depends on the 
gluon momentum as well as on the vector index of the gluon and the magnetic quantum 
number of the P-state. 

In the conclusions of sect. 5 we briefly discuss how an extension of our perturbative 
results to a wider range of applicability can be carried out. Applications to P-state 
annihilation decays need some further effort and will be presented elsewhere. 

2 Coupled Bethe-Salpeter equations 

We consider a mesonic bound state \K > with four momentum . We are interested 
both in the quark- antiquark amplitude 

-iKX . 

^ QQ (x 1 ,x 2 ;K) =< 0\T^( Xl )^(x 2 ))\K >= J e^^fo K)d 4 p, (2.1) 

as well as in the quark-glue-antiquark component 

^ qa ^q(xuX 3 ,x 2 ;K) = < Q\T{^{x 1 )At x {x^{x 2 ))\K > 

-iKX 

= -j^s J e- l{pr+qx) ^QA»Q(p,q-,K)d 4 pd\ (2.2) 

Spin- and color indices are suppressed, explicitly the color singlet states are understood 
as # = 0-/VN)M a , and ipA^ifj = ^2/{N 2 - l)^ a2 A» a {\ a a2ai /2)^ ai , with N = 3 the 
number of colors. There is no need of specifying the gauge at this stage. 

In the two particle QQ system, with masses mi,m 2 , we used the usual CM and 
relative variables in position and momentum space: 



X = rfiXi + i] 2 x 2 , r = xi - x 2 ; x 1 = X + rj 2 r, x 2 = X - r} x r, 
K = pi+p 2 , p = T] 2 pi - 7] X p 2 ; p 1 = r] 1 K + p, p 2 = r] 2 K - p. (2.3) 

Here rjj = mj/{m\ + m 2 ), the total and reduced mass are denoted by M = mi + m 2 
and n = m\m 2 j[m\ + m 2 ). For the three particle system QAQ it is convenient to use 
the same CM and relative coordinates X, r for Q and Q as above. We denote by x the 
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distance of the gluon to the CM of the QQ, the momentum of the gluon is denoted by 
q. The coordinates then read 

x = x 3 - 771X1 - r] 2 x 2 ; x 3 = X + x, 

K = px+p 2 + q ; p! = r)i(K - q) + p, p 2 = r) 2 (K - q) - p. (2.4) 

The derivation of the bound state equations for iPqq(p; K) and iPqa»q(p,Q] K) is a bit 
technical. One may skip the rest of this section, the resulting lowest order graphs in Fig. 
6 and Fig. 7 of the next section are easily understood anyhow. For a systematic and 
general approach we consider the Green functions 

G Q q jQQ (x 1: x 2 ,x' 1: x / 2 ) =< 0\T([^(xMx 2 )Mx[)^(x' 2 )])\0 > (2.5) 




- i(pr+p ' r ' +K{X+X,)) G QQ , QQ (p, p'; K)d'pd 4 p'd 4 K, 



and 

(x 1: x 3: x 2: x' 1: x' 2 ) =< 0|T([V(xi)^(x 3 )^(x 2 )][^(x / 1 )^(4)])|0 > (2.6) 

Further Green functions, which will be needed below, are defined in an analogous way 
and immediately recognized from their symbols. Wave functions and Green functions 
in momentum space are denoted by open circles with the appropriate legs attached to 
them. The quark gluon vertex function is defined in the usual way as the sum of all 
one-particle irreducible graphs contributing to the QAQ-vertex, i.e. the graphs which 
cannot be separated into two different pieces by cutting only one internal line. The 
external legs are not included. 

With the help of the vertex function the Green function Gq A ^qqq can then be 
decomposed as in Fig. 1, where internal lines, here as well as in the other graphs of this 
section, denote full propagators. 




Figure 1: The decomposition of Gq A ^qqq. 

We define the totally connected Green function G c ^k qq, denoted with the filled 
circle in Fig. 1, as the sum of the graphs which cannot be divided into two separate 
pieces by cutting one internal Q (or Q) line plus the left external Q (or Q) line. 

In order to separate the different components in the bound state we have to intro- 
duce four irreducible kernels Iqq,qqJqq,qa»qJqa»q,qqJqavq,qa»q- They are defined 
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as the sums of all connected graphs contributing to the corresponding Green functions 
Gqq qq, • • • , respectively, without external legs, and will be denoted by rectangles with 
the appropriate lines attached. The irreducible kernels have the following properties: 

a) The graphs cannot be separated into two pieces by cutting two internal lines of the 

type QQ, or by cutting three internal lines of the type QAQ. 

The kernels Iqq,qa»q an d Iqa»q,qa v q have a QAQ state on the right hand side, and 
there are three further rules for them: 

b) The graph cannot be separated into two pieces by cutting an internal Q (or Q) plus 

the right external Q (or Q) line. 

c) The graph cannot be separated into two pieces by cutting two internal QA (or AQ) 

lines plus the right external Q (or Q) line. 

d) The graph cannot be separated into two pieces by cutting two internal QQ lines plus 

the right external A line. 

The prescriptions b), c), d) only refer to lines on the right hand side. If the 
corresponding situation appears on the left hand side, the graph is included in the 
definition of I qa ^q,qq or I qa »q,qa»q- 

The above definitions have been chosen in such a way that one is able to extract 
the desired wave functions while avoiding double counting in the decompositions below. 
Note that we do not make use of the usual two quark irreducible kernel / consisting 
of all connected graphs for Gqq qq (without external legs) which cannot be separated 
into two pieces by a QQ cut. We need a more detailed classification of kernels for our 
purpose. 

Consider now the sum of all graphs which contribute to Gqq qq. They may be 
collected in different classes using the above definitions. This is shown in Fig. 2. 




Figure 2: The equation for G Q q <Q q. 
The corresponding decomposition of Gq A ^q qq is shown in Fig. 3. 
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Figure 3: The equation for G qa ^q Q q. 

The correctness of the decompositions in Fig. 2 and Fig. 3 can be demonstrated in 
any order of perturbation theory by showing that there is a one to one correspondence 
between graphs. Clearly each graph appearing on the rhs also belongs to the lhs. The 
proof that each graph on the lhs appears exactly once on the rhs is given in the ap- 
pendix. There we also show some examples for graphs which are forbidden by our above 
definitions. 

The intermediate state \K >< K\ leads to a pole in the Green functions at K° = 
±(yj ' M\ + K 2 — it) on both sides of Fig. 2 and Fig. 3. Extracting the residue and 
dropping the common factor ~ ^qq we arr i ve a ^ ^ wo coupled Bethe-Salpeter equations 
for the wave functions ipQQ and iPqa^q which are displayed in Fig. 4 and Fig. 5. 




Figure 4: The equation for ipQQ- 
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Figure 5: The equation for i/jqa^q- 

3 Perturbative solution 

Up to now the treatment was general and there was no need to specify the Lagrangian or 
the gauge in detail. In this section we will specialize to non relativistic QCD (NRQCD 
up to order 1/m) with the Lagrangian j2] 



Here ip is the Pauli spinor for the quark field (annihilation operator), and \ f° r the 
antiquark field (creation operator). Our conventions are _D M = d^ + igA^, F^ v = d^A v — 
d v A» + ig[A», A% E n = -F 0n , B n = - e n i k F jk /2. Since we are working in lowest order 
only, there is no need of introducing matching coefficients in the Lagrangian (3.1). We 
go to the rest frame of the bound state where = (— \E\, 0). Note that the rest masses 
are dropped in NRQCD, we will also suppress the K in our notation from now on. The 
perturbation theory is performed in radiation gauge VA = 0, the gauge best suited for 
treating bound states. From now on we distinguish between Coulomb gluons A denoted 
by dashed lines, and transverse gluons A denoted by curly lines. 

We first solve the bound state equation in Fig. 4. The leading contribution to the 
kernel Iqq^qq which appears in the first term on the rhs is given by the exchange of one 
Coulomb gluon between Q and Q which leads to a static kernel: 



»(27r)V(A^/2)(A^/2) 

inn nn — ; rr^ • {^■^) 



l QQ,QQ - | p _ p ,|2 

The corresponding contribution with a transverse gluon contains additional factors of 
p/ra at the vertices. It is therefore suppressed if we consider only soft external momenta. 
Internal momenta are also essentially soft due to the non relativistic wave function given 
below. The last three contributions in Fig. 4 are also of higher order. In lowest order 
Fig. 4 therefore gives the familiar non relativistic Bethe Salpeter equation in ladder 
approximation and radiation gauge. Because the kernel is static, the only dependence 
on the relative energy variable p° on the rhs is in the propagators. One can therefore 
integrate over p° on both sides and derive the Schrodinger equation for the momentum 
space wave function ^qq(p) = (27r) -5 / 2 / ^qqIPiP )^ m the usual way. After this has 
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been solved, the dependence on the relative energy can be recovered by going back to 
the original equation. The well known solution then reads 



i(2n) 


3/2 (\En, 


+ p 2 1 /2m 1 + p 2 2 /2m 2 )^nim(p) 


(p° - 77x | 


E n , 


- p 2 /2m 1 + ie)(p° + r) 2 \E n \ 


+ p 2 2 /2m 2 - ie) 



( I _\ i 0\ _ v\^n) Vl-^ral T I'll ""'n T \>2l A ""'2-)H J nlm\V) / Q q n 

{VQQ)nlm{P, V ) ~ / n „ |pi |_ n 2 rtra , 77 I T „2 /o„ " ^ 



The spin wave functions factorize. i? n = — (CVa: s ) 2 /i/ (2n 2 ), with CV = (iV 2 — 1)/ (2iV) = 
4/3 the Casimir operator in the fundamental representation, is the energy eigenvalue, 
and ipni m (p) = VVtf(|p|)^m(p) the normalized Schrodinger wave function in momentum 
space. The normalization has been chosen such that / ip(p,p )dp = (27r) 5 / 2 -^(p), which 
means that ipqQ in (2.1) becomes the correctly normalized Schrodinger wave function 

if x\ = x®- The Bohr radius is a = 1/ ^j2^\E~y\. For the lowest states one has (it is 
convenient to write them in terms of the energy E n instead of the Bohr radius) 



V2(2//|.Ei 



^5/4 



<MIp|) = 

r n n 8y^(2 / .|E 2 |) 5 / 4 (|p| 2 -2 / .|E 2 |) 

W|P|) = x— |o I .in ixa ' ( 3 - 4 ) 



<M|p| 



7(p 2 + 2/^1^1 

^|^ 2 |) 5/4 ( 

v^dpl 2 + 2n\E 2 \f 
16v / 2(2 At |£7 2 |) 7 / 4 |p| 



37r(|p| 2 + 2^1^1)3 



Having solved the bound state equation in Fig. 4 in lowest order, one can insert the solu- 
tion, together with the leading approximations for the propagators, the vertex functions, 
and the kernels, into the rhs of Fig. 5, dropping the higher order contributions of the 
second line. The power counting in the coupling constant g requires some care, because 
there are not only the explicit factors g in the vertices, but also factors p/m ~ g 2 in the 
coupling of transverse gluons to quarks, as well as momentum factors in the three gluon 
vertex. Taking this into account, the leading contribution to tPqa°q is then given by Fig. 
6, the leading contribution to ^qaq by Fig. 7. 

p' 

q —-( j= = q C j= + q 
p — P P 

2 2 2 

(a J 




Figure 6: The leading contributions to iPqaoq- 

In Fig. 7 c) the kernel Iqaq,qq has been replaced by its lowest order contribution, 
where the transverse gluon couples to two Coulomb gluons, while the latter couple to 
the quark and antiquark, respectively. For the reasons mentioned above, this graph, 
although of higher order in the coupling constant, contributes to the same order as the 
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Figure 7: The leading contributions to V'qaq- 

tree graphs in Fig. 7 b), because the direct coupling of the transverse gluon to the 
quark or antiquark in 7 b) is suppressed by a factor p/m ~ g 2 . For the Coulomb gluon 
the corresponding graph does not exist, because there is no coupling of three Coulomb 
gluons. There is no non abelian contribution in leading order in Fig. 6. 

We can now insert the QQ wave function (3.3) into the rhs of Fig. 6 and Fig. 7, 
and calculate the wave functions on the lhs. We start with the tree graphs Fig. 6 a) 
and Fig. 7 b). They can be simply written down because there is no integration over 
internal momenta. We are interested in the equal time wave function, therefore we first 
integrate over p° (or, more simply, over p 2 m a i),bi) and over p\ in a 2 ),b 2 )) with the 
residue theorem. There are three denominators, one arising from the external Q or Q 
propagator and two from the QQ wave function. For ai),&i) e.g. one gets 



/ 



dpi 



(p° 2 + q° + \E n \ + p?/2mi - ie)(p° 2 + \E n \ + p'?/2m 1 - ie)(p° 2 - p|/2m 2 + ie) 

= ^ f35 ) 

(q° + \E n \ + p?/2mi + p|/2m 2 - ie)(\E n \ + p' 2 1 /2m 1 + p|/2m 2 ) ' V ' ; 

Next we have to perform the integration over q°. For the Coulomb gluon the propagator 
is independent of q°, one thus cannot simply perform the g°-integration, because the 
denominator in (3.5) contains only one power of q°. We therefore choose an infinitesimal 
relative time t of the gluon with respect to the quark and antiquark. The Fourier 
transform then becomes 



/ 



e -> q °t dq o 



2-KiQ(—t) for t infinitesimal. (3.6) 



(q° + \E n \ + p?/2mi + pi/2m 2 - ie) 



We will comment this result below. Finally we add the contribution of Fig. 6 a 2 ) and 
divide by (27r) 5 which gives the correctly normalized equal time wave function ip(p,q) 
in momentum space. This results in 

<Wp> q) = 9 ^,l\J (^(-p 2 ) - ^(pO) • (3-7) 

For a transverse gluon the propagator provides enough powers of q° in the denominator 
to make the equal time wave function well defined. The result from the graphs Fig. 7 
bi), b 2 ) becomes 
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(^QAQ) (6) (P,q) = 7 9 ~^- T>< 

Q Q (27r)3/22|q|(|q| + |E n |+pf/2m 1 + p|/2m 2 ) 

( [(Pi)r + 'qx Si]^(-p 2 ) [(p 2 )r + x s 2 ]^( Pl ) \ 
V mi rn 2 ) ' 

Here p T = p — (pq)q/q 2 denotes the transverse component of p. 

We finally come to the non abelian graph Fig. 7 c). This involves an integration 
over the internal four momentum p' . Performing first the integrations over p'°,p°, and 
q° gives the equal time wave function 



0W (c) (p,q) = 7 aNy/u * -x 

Q Q (27r)9/2 2 |q| (|q| + \E n \ + pf/2m 1 + p|/2m 2 ) 



/ 



(p - p')MpW 



p - p' - ^iq| 2 |p - p' + ^ql 2 



We come to the physical interpretation of the above results. First we have to recall that 
there are three different scales in NRQCD: Hard (~ m) ^> soft (~ a s m) ^> ultrasoft 
(~ a 2 s m). Quark masses are hard, quark momenta soft, and quark kinetic energies as 
well as binding energies are ultrasoft. Since we used NRQCD, the momenta pi,P2,q in 
the QAQ wave functions (3.7) - (3.9) must also be soft. In fact these expressions are 
dominated by soft momenta. 

The graphs in Fig. 6 describe the Coulomb field moving along with the quark and 
antiquark. It is easy to perform the Fourier transform of (3.7) with respect to p and q. 
It becomes 

One recognizes the product of the QQ-wave function with the Coulomb potential gen- 
erated by the Q and the Q source. (In QED one has, of course, to replace gy/Cp by 
e.) The appearance of the step function 0(— t) is due to the time ordering prescription 
in the wave function (2.2). For t > the gluon field A°(xs) has to stand to the left. 
One can express it in terms of the quark charge densities, and the latter annihilate the 
vacuum in lowest order. 

The graphs in Fig. 7 b) describe the fields which arise because the quarks are moving. 
For large |q| the spin independent contributions of the quark and antiquark are those 
of the Coulomb contributions (3.7) for t < 0, multiplied with (1 /2)(p j) T /rrij. With 
the exception of the factor 1/2 in front, this is exactly what one would expect from 
classical electrodynamics. The appearance of this factor can be understood as follows. 
In classical radiation problems one uses the retarded Green function. In field theory, like 
in our definition of the wave function (2.2), on uses, however, the time ordered product. 
If we had used the retarded product instead, the denominator of the transverse gluon 
propagator would have been changed from q 2 + ie to (q° + ie) 2 — |q| 2 , such that both 
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poles in the complex g°-plane would lie in the lower half plane. In this case one would 
simply close the contour around the remaining g°-pole from the quark propagator in the 
upper half plane and would end up with the replacement 



2|q|(|q| + |£ n |+p 2 /2m 1 + p 2 /2m 2 ) ' | q |2 _ (|£ n | + v \j2m x + p|/2m 2 - tef ' 

(3.11) 

removing the factor 1/2 for large |q| in (3.8). The use of the advanced propagator would, 
on the other hand, give a vanishing result. This shows that one has to define the equal 
time wave function very carefully. 

In the following we will stick to our definition with the time ordered product. If 
one defines the equal time wave function for A as the average of the expressions for 
positive and negative infinitesimal t, the replacement 0(— t) =^1/2 leads to the correct 
relation between A and A for large |q|. For smaller |q| the result for A is modified, 
the denominator is less singular for |q| — > 0. This is, of course, due to the fact that the 
quarks are not free, there is an additional long range radiation field. Fig. 7 c) finally 
arises from the glue moving along with the quarks, it has no simple classical analog. 

The expressions for the transverse gluon can be simplified. From the denominators 
in (3.8) and (3.9) it is obvious that these expressions are dominated by ultrasoft q. If we 
thus restrict to ultrasoft q and remember that pi and p 2 are soft, we may approximate 
(3.8) by using pi = p - 771 q « p, -p 2 = p + r] 2 q « p, and pf/2m 1 + p 2 /2m 2 = 
p 2 /2/i + q 2 /2M « p 2 /2/i, and get 

0W (6) (p,q) = 3VC~ f pMp) 

Q Q (27r)3/2 2/ i|q| |q| + KI+ P 72//) 



We dropped the spin dependent term because it is suppressed for ultrasoft q. 

We next discuss the non abelian contribution (3.9). For ultrasoft q one may again 
simplify the kinetic energies in the denominator, furthermore one can drop the q in 
the integrand. The integral can then be easily performed for a given ip(p'), e.g. by 
introducing a Feynman parameter. We give the results for the three lowest states. 

(^QAQ)lOo(P.q) = , n 7Z. 77, , ~ ;o s, Z ~ r= ,, ^00, (3-13) 



(27r)3|q|| 
g 3 Ny/i 


;|q| + |£; 1 |+pV2//)(p 2 + 2 / /| J E; 1 |) 2 
^(2 /U |E 2 |) 3 / 4 (p 2 -3-2 / u|E 2 |)p T 


(27r) 3 |q|| 




+ 


E 2 \ 


+ p 2 /2 /U )(p 2 + 2 / i|E 2 |)3 



(^AQ) 2 oo(p,q) = \:::~ r, :r: v :,. r. r' 1 ? 1 :::'.., ^- ••• | o 

The result for the P-state consists of two parts: 



(4>, 



QAQ)21m 



vP,q) 



Ug 3 N^/C^(2fi\E 2 \) 5 / 4 



x 



(27r)3 v ^|q|(|q| + | J E 2 |+p2/2 /U ) 

PT|p|yim(p) _ 1 (e fc ) r Yj m (e fc ) 
p 2 + 2^1^1)3 4 (p 2 + 2/i|E 2 |) 2 



(3.15) 
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where denote the three Cartesian unit vectors and (e^r their transverse parts. 

For the ground state, as well as for the first term of the P-state, we may express the 
result directly in terms of the wave function and compare it with the result of Fig. 7 b) in 
(3.12). The non abelian contribution (3.13) for the ground state is N/ (2Cp) = 9/8 times 
the result in (3.12), while for the P-state the first contribution in (3.15) is N/Cp = 9/4 
times the corresponding result in (3.12). Thus numerically the non abelian loop graph 
is more important than the tree graphs. 

In momentum space the leading order of the transverse electrical field = — d°A 
and the magnetic field B = V x A are obtained by replacing V => iq and d° —iqo => 
— i\q\ and applying these factors in the expressions for A. The last replacement is due 
to the fact that the g°-integration was performed by closing the contour around the pole 
at q° = |q| — ie. The longitudinal electric field = — VA° can be obtained from (3.7) 
by multiplying with — iq, or directly in position space from (3.10). 



4 The QAQ wave function of the P-state at r = 

A quantity of particular interest is the amplitude for the special case of quarks with 
identical flavor at zero distance in the QAQ component of the P-state. This component, 
although suppressed, contributes to annihilation decays to the same order as the QQ 
annihilation in the leading QQ component. The reason is that this annihilation is also 
suppressed, because the wave function at the origin vanishes for P-states. On the other 
hand the QQ in the QAQ wave function can be in an 5-state. In one loop order one 
finds that logarithmic divergences are canceled by adding up both contributions. This 
has been discussed in detail by Bodwin, Braaten and Lepage PQ. Let us therefore consider 
the wave function for quarks at zero distance 

^QAQ(q)m = (27r)- 3/2 J 4> QAQ (p, q) 2 lmd 3 p. (4.1) 

The calculation is simple if we restrict to the dominating region where q is ultrasoft. 
From (3.12), (3.15) we find 



7 , v ig^{m\E 2 \fl\e k ) T Y lm {e k ) 

WQAQ\SL)m - = ; 3 X 

3\/67r 3 / 2 |q| (ym(|q| + \E 2 \) + ^m\E 2 



2N 



3^m(|q| + \E 2 \) + \Jm\E 2 \ - — ^m\E 2 \j . (4.2) 

The first two terms in the second line originate from the graphs in Fig. 7 b) and the 
resulting contribution (3.12), the third one is a sum of (n/C f ) {3^m(|q| + \E 2 \) + 

^m\E 2 \} from the first part of (3.15) and of -(3N/C F ){ v /m(|q| + \E 2 \) + ^fm\E 2 \} 
from the second part. 

Our result (4.2) shows a remarkable feature. Recalling that 2N/Cf = 9/2, one finds 
that the non abelian term dominates the abelian one for small |q| and makes the second 
line of (4.2) negative. For |q| = (13 /36) I-E2I the bracket vanishes, for larger |q| the 
positive abelian term dominates. 
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The wave function (4.2) determines the amplitude for annihilation of the two quarks 
in the color octet state into gluons which subsequently hadronize. As one might have 
expected, (4.2) depends on the gluon momentum q, as well as on the relative orientation 
between the gluon and the angular momentum of the P-state. It is, however, independent 
of the quark spins. The QQ wave function has entered in a complicated way. It is not 
possible to express the result simply by the derivative -R'(O) of the radial wave function 
at the origin. 

5 Conclusions 

In the foregoing we derived some results which may appear rather unexpected at first 
sight. One of them is the discontinuity of the wave function in the relative time of 
the Coulomb gluon with respect to the quarks. The other one is the relation between 
transverse and Coulomb gluon. Both results could, however, be understood, the point is 
the use of different Green functions (retarded versus time ordered) in classical radiation 
theory and in field theory. 

One may wonder whether some of these features could be artifacts of the gauge 
or of the perturbative treatment. We are confident that this is not the case. One 
could, of course, choose manifestly gauge invariant wave functions by connecting QQ 
or QA^Q, respectively, with path ordered strings. The QQ wave function would now 
contain already some glue in it due to the presence of the string. This glue along the 
straight string is, however, not an additional degree of freedom. Everything is fixed by 
the positions of the quark and the antiquark. The essential feature of a QA^Q wave 
function is not the presence of some glue, but the presence of an additional degree of 
freedom carried by the gluon. 

The strings would not contribute in lowest order perturbation theory. A consideration 
of higher order perturbative contributions within the approach presented here, whether 
with or without strings, is certainly possible. It will not alter the essentials of our results. 
Our method could also be extended to obtain wave functions for multi gluon admixtures 
which arise in higher orders. 

A systematic treatment of non perturbative effects is beyond the present investiga- 
tion. We emphasize, however, that our whole approach can be immediately general- 
ized to the case where we replace the kernel (3.2) by some other static kernel, say the 
Fourier transform of a phenomenological static potential involving confinement. This 
will, in turn, modify the QQ Schrodinger wave function and the subsequent results. For 
phenomenological applications such an approach appears quite legitimate, it has been 
familiar and successful in other fields of quarkonium physics. 

The application of our results to production or decays of heavy quarkonia goes be- 
yond the present investigation. In the framework of NRQCD, hard annihilation processes 
are described in lowest order by four fermion terms of dimension six. They were not in- 
cluded in our Lagrangian because they are of order 1/m 2 and thus not important for 
the calculation of the wave function. While the hard subprocess of «S-wave annihilation 
is simple and well understood, the complicated dependence of our wave function on the 
gluon momentum and the magnetic quantum number of the P-state require some addi- 
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tional work on this topic. While we know that logarithmic divergences which show up in 
higher orders cancel exactly against corresponding divergences in the P-wave annihila- 
tion of the QQ component, there is no information about the finite part up to now. This 
information can be obtained from our explicit wave function. Due to the compensation 
between the two contributions in (4.2) and the change in sign between smaller and larger 
|q|, suppressions compared to naive expectations may be expected. The results will be 
presented elsewhere. 

Acknowledgement: I thank D. Melikhov and N. Brambilla for reading the manuscript 
and for valuable suggestions. 

References 

[1] G. T. Bodwin, E. Braaten, G. P. Lepage, Phys. Rev. D 46, R1914 (1992); 
with T. C. Yuan, Phys. Rev. D 46, R3703 (1992); 
Phys. Rev. D 51, 1125 (1995), Phys. Rev. D 55, 5853(E) (1997). 

[2] W. E. Caswell, G. P. Lepage, Phys. Lett. B 167, 437 (1986). 

A Decomposition of Green functions 

In this appendix we prove the correctness of the decompositions in Fig. 2 (the proof for 
Fig. 3 is analogous) by showing that every graph on the lhs appears exactly once on the 
rhs. Consider an arbitrary graph contributing to the lhs. 

1) If the graph is disconnected it uniquely belongs to the first term on the rhs of the 
first line of Fig. 2. 

Consider a connected graph next. 

2) If the graph cannot be separated by a QQ or a QAQ internal cut, it belongs to 
InnnQ plus external lines. It then belongs to the second term on the rhs, where the 
disconnected part of the Green function appears. 

Let's next take the case where the graph can be cut by an internal QQ or QAQ cut. 
The propagators are always taken as part of the right hand graph. We first construct a 
unique preliminary cut C, p as follows: Follow the quark line from the left until one comes 
to an internal propagator where the graph can be cut in the above way. If there are 
several possibilities to do so, let C p be the cut where the antiquark line is cut as far to 
the left as possible. 

3) If C p is a QQ cut we take it as the final cut C. In this case the graph belongs 
to the second term on the rhs of Fig. 2 with a connected contribution to the Green 
function. 

If C p is a QAQ cut there are two alternatives: 

4) It may be possible that one can transform the cut by keeping it's position at the 
Q (or Q) line, while moving the position of the cut at the Q (or Q) to the right, such 
that one obtains a QQ cut. This happens if one can shift the cut to the right by moving 
it over a vertex part, such that it does no longer cut the gluon line. If this procedure is 
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possible at all, it is unique. In this case we take this QQ cut as the final C. Note that 
the graph left to the cut is a contribution to Iqq^qq; it is not allowed to cut the graph 
by the original C p , because it would cut an external Q or Q line now. Again we get a 
contribution to the second term on the rhs of fig. 2. 

5) If the above manipulation is not possible, we take the original QAQ cut C p as the 
final C. In this case the graph on the lhs is a part of Iqq qaQi while the one on the rhs 
is a part of the totally connected Green function G C qa%qq- Contributions to the first 
two terms on the rhs of Fig. 1 cannot appear by construction because our prescription 
would have moved the corresponding cut to the right. 

We have thus shown that every graph on the lhs also appears on the rhs. In order 
to proof that it belongs exactly to one term on the rhs we will next show that it is 
not possible to move the cut C to the right, producing thereby another admissible cut 
C r , which would lead to a different decomposition. (By construction it is certainly not 
possible to move the cut to the left). 

Consider first the case that C is a QQ cut. If we move it to the right, crossing a 
vertex both at the Q line as well as at the Q line (irrespective whether we arrive at an 
admissible cut) , the remaining graph on the left could be cut by the original C and thus 
would not belong to an irreducible kernel. If we move the cut only over one vertex, such 
that it becomes a QAQ cut, the graph on the lhs would be inacceptable because of rule 
b) in the definition of the kernels in sect. 2. 

Consider next the case that C is a QAQ cut. If we move it to the right by crossing a 
vertex both at the Q line as well as at the Q line (irrespective whether we arrive at an 
admissable cut), the remaining graph on the left could be cut by the original C and thus 
would not belong to an irreducible kernel. Rule d) of sect. 2 applies in the case that the 
gluon line which was cut by C is also cut by C r and thus would become an external line. 

If we move the cut only along the Q (or the Q) line while keeping the position at the 
Q (or the Q) line we may either obtain a QQ cut or a QAQ cut (or an inadmissible cut). 
In the first case the kernel to the left is forbidden by rule b), in the second case by c). 

The somewhat complicated definition of the kernels Iqq qaqi Iqaq,qaq with respect 
to external lines at the rhs was necessary to guarantee the uniqueness of the decompo- 
sition. In Fig. 8 we show some graphs which do not belong to Iqq qq or to Iqq : qaq, 
together with the cuts and the rule of sect. 2 which prevents this. 




a) b) c) d) 

Figure 8: Some reducible kernels. 
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